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Abstract. We develop an effective low-frequency theory of the electromagnetic field in equilibrium with thermal 
objects. The aim is to compute thermal magnetic noise spectra close to metallic microstructures. We focus on the 
limit where the material response is characterized by the electric conductivity. At the boundary between empty space 
and metallic microstructures, a large jump occurs in the dielectric function which leads to a partial screening of low- 
frequency magnetic fields generated by thermal current fluctuations. We resolve a discrepancy between two approaches 
used in the past to compute magnetic field noise spectra close to microstructured materials. 

PACS. 05.40. -a Fluctuation phenomena and noise - 03.75. Be Atom optics 



1 Introduction 

In the context of atom chips, low-frequency thermal magnetic 
noise has recently emerged as one crucial element that limits 
the lifetime of miniaturized atom traps. Recent experiments 
I1I2I3I have confirmed the basic features of theoretical pre- 
dictions for spin flip processes induced by magnetic near field 
fluctuations. A current trend is to extract other physical mecha- 
nisms that lead to loss by subtracting the near-field induced loss 
rate. One particularly interesting mechanism is the lowering of 
the trap depth due to atom-surface interactions j3). Accurate 
calculations of magnetic near field noise are clearly needed for 
this purpose. Magnetic fluctuations are also relevant in other 
contexts, for example in biophysics where they impose ulti- 
mate limits on the sensitivity of SQUID detectors |4|, and in 
magnetic resonance force microscopy, a near-field variant of 
magnetic resonance imaging 15 161 . 

Typically, one is interested in field frequencies hco -C ksT 
where the noise is dominantly classical. The border to the quan- 
tum regime can be reached with magnetically trapped atoms, 
either by cooling the microtrap components and/or applying 
strong static magnetic fields that push up the relevant frequency 
range (given by the Larmor frequency). Even at the highest 
frequencies conceivable with state-of-the-art atom chip struc- 
tures (in the GHz range), the (vacuum) field wavelength A is 
much larger than the characteristic distances, so that the qua- 
sistatic approximation applies outside the structures. This leads 
to a peculiar situation to describe the field fluctuations. Indeed, 
one cannot apply the standard procedure and attribute thermal 
or quantum fluctuations to the normal mode amplitudes of the 
field, because there are no nontrivial solutions to the homoge- 
neous field equations (i.e., eigenmodes) in the quasistatic limit. 
Near field noise is actually dominated by the fluctuations of 



its sources (currents, magnetic moments) whose spectral mode 
density depends on material or atomic constants I4I7I8I9I10I . 
As a result, the near field noise spectrum differs markedly from 
the celebrated blackbody radiation law 1 411 11121 . 

Roughly, two approaches can be identified to compute mag- 
netic noise close to micro- and nanostructures. The first one can 
be traced back to the fluctuation electrodynamics put forward 
by Rytov and co-workers 1131 in the 1950's. Based on a sta- 
tistical thermodynamics argument (the fluctuation-dissipation 
- FD - theorem 1141 ). random charge and current fluctuations 
are associated to a dissipative material structure. Their radia- 
tion is incoherently summed to give the total noise strength of 
the field. In a planar geometry, the radiated field and the re- 
quired averaging can be calculated analytically. Results along 
this line have been computed and experimentally verified for 
planar metallic layers by Varpula and Poutanen in 1984 |4|. 
Sidles and co-workers give an extensive discussion with appli- 
cations for magnetic resonance microscopy and quantum com- 
puting |6|. 

An alternative approach uses the FD theorem for the (mag- 
netic) field itself and has been popularized in a series of papers 
by Agarwal in 1975 1 15 1. The advantage is that the incoherent 
averaging is avoided; the FD theorem reduces the calculation 
to the radiation of a single dipole source (Green function), lo- 
cated at the observation point. This method has been applied, 
in the context of atomic microtraps, by the present author and 
co-workers 1 16 1 and Rekdal and co-workers 11711 81 . Both ap- 
proaches have been shown to be equivalent under fairly general 
conditions, thanks to an identity that implements the FD theo- 
rem in electromagnetism 1 191201 . 

In the context of integrated atom optics, incoherent sum- 
mation over fields has been put forward by S. Potting and the 



2 



Carsten Henkel: Magnetostatic field noise near metallic surfaces 



present author as a versatile tool to handle arbitrary nanostruc- 
tures 1211 . It just remains to perform a certain spatial integral 
over the volume filled with electrically conducting material. 
This yields the correct scaling of the noise spectrum with the 
atom chip geometry, provided the skin depth is long enough. 
Based on this approach, for example, the Vuletic group could 
fit reasonably well spin flip loss rates close to rectangular wires 
1 3 1. (See also Refs. 1 18 22 1 for a re-analysis and discussion.) A 
closer comparison shows, however, that the theoretical results 
are off by numerical factors between two and three compared to 
the noise spectrum predicted by the FD theorem |21 1. This dis- 
crepancy is the motivation for the present paper. We point out 
an error in the 'incoherent summation' approach that is linked 
to the particular boundary conditions for the electromagnetic 
field at the surface of a good conductor. We derive approximate 
boundary conditions that apply to any geometry in the low- 
frequency range relevant for atom chips. In the planar case, we 
show that they lead to an accurate agreement with the 'Green 
function' approach in the limit that the vacuum wavelength is 
the largest length scale. The only point missing in the theory 
is the blackbody noise level that prevails at large distances, but 
this one is in most situations impossible to detect anyway. 

The parameter regime we focus on in this paper is illus- 
trated in Figure^ We shall call 'quasi-static' the regime where 
the skin depth 5 inside the material is larger than any other 
geometrical scale (denoted a in FigQ. We focus on metal- 
lic materials and use the definition S = in 
terms of the (DC) conductivity a. Our theory aims at cover- 
ing both the quasi-static regime and a skin depth comparable 
to a. The temperature T defines another frequency scale be- 
low which the field fluctuations behave classically. This is ac- 
tually not a limitation as long as we assume thermal equilib- 
rium. The theory is extended into the quantum regime with 
the replacement fc^T i— > ^hui coth(huj /2ksT) and assum- 
ing symmetrized noise correlation spectra. At frequencies in 
the visible and ultraviolet range, however, the dielectric func- 
tion of the material becomes complicated, and more parameters 
(transverse optical phonon frequency, plasma frequency, tabu- 
lated data . . . ) are needed for an accurate modelling. 

The present paper thus aims at clarifying the validity of 
the 'incoherent summation' approach and at extending it into 
the regime of a short skin depth using the appropriate bound- 
ary conditions. We also argue that from a practical point of 
view, the FD approach appears simpler because it is sufficient 
to compute the radiation from a single point source, while for 
'incoherent summation', many sources (anywhere inside the 
spatial domains filled with absorbing material) have to treated. 

In the following sections, we start by writing down the ba- 
sic equations for the magnetic field and explain the relevant 
parameter regimes (Sec|2j. The boundary condition at the sur- 
face of a good conductor is derived. In section|3] we review the 
incoherent summation technique and show for the special case 
of a metallic half-space that with the correct boundary condi- 
tion, one gets a magnetic noise spectrum in agreement with the 
FD approach. We show that in the limit of a short skin depth, 
the transmission of the magnetic field out of the metal becomes 
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Fig. 1. Characteristic frequencies involved in near field electromag- 
netic noise. The conductivity of silver, <7Ag is used as a convenient 
scaling parameter. A typical geometrical feature size is denoted by 
a. Three characteristic frequencies separate different regimes as illus- 
trated by the vertical lines. The formulas at the frequency axis give 
explicit values and their scaling with the relevant parameters. 
For example, on a scale a = 10 /im, the skin effect is irrelevant 
for the field propagation near a silver structure at frequencies below 
0.1 MHz x 10 4 = 10 3 MHz; the quasi-static approximation applies in 
this regime. At frequencies approaching 2.5 10 17 Hz, silver becomes 
transparent, and its permittivity does no longer involve a purely real 
conductivity. 



much less efficient. In subsection !3.4l we give a qualitative ex- 
planation of the power laws in the distance-dependence of the 
noise spectrum and review results obtained for a thin metal- 
lic layer. In sectionf4] we formulate the equations to be solved 
within the FD approach when the quasistatic approximation is 
made in the spatial domains filled with vacuum. The formula- 
tion applies to an arbitrary geometry and is then specialized to 
a metallic half-space. In the latter case, we demonstrate agree- 
ment with the more complex, fully retarded FD approach in the 
long vacuum wavelength limit. 



2 Boundary conditions at low frequency 

We want to solve the Maxwell equations in a non-magnetic 
medium (/i(x) = /i U ): 

V-(oE+j) = 0, VB = 0, (1) 
V x E = iwB, V x B = fig (crE + j + V x M) .(2) 

where cr(x) is the metal conductivity. In the vacuum regions, 
cr(x) = — isqu) which will be assumed much smaller in mag- 
nitude than the metal conductivity, denoted a. Two kinds of 
sources appear here that apply to the methods of magnetic 
noise calculations mentioned in the introduction. They are rep- 
resented by the externally prescribed terms j(x) and M(x). In 
the incoherent summation technique, the current density j (x) 
is localized inside the metal and represents thermal fluctua- 
tions. It is a random quantity with correlation function given 
by Ea.( ll4> . For the calculation of the magnetic Green function, 
the magnetization M(x) can be identified with the magnetic 
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moment of an atom localized in vacuum outside the metal, see 
after Ea A37l . Depending on the method of calculation, only 
one or the other source term is actually nonzero. 

Taking the curl of l^), we find the following wave equation 
for the magnetic field 



v 2 b 



? 2 B 



-Mo [V x j(x) 
+ (Vcr) x E] 



V x V x M 



(3) 



Outside the metal, g(x) = uj/c. Inside the metal, q(x) is com- 
plex and related to the skin depth 8 by q 2 = 2i/5 2 = i[io<ju. 
We could also introduce a spatial diffusion coefficient D = 
Hoa by going back to time-dependent equations. The source 
terms of the magnetic wave equation Q are related to the cur- 
rent density, the magnetization density and to a current parallel 
to the surface induced by the electric field. It turns out that this 
last current leads to a jump in the derivative of the magnetic 
field across the interface. The solutions to the wave equation 
can be sought in terms of a Green function, assuming local- 
ized sources. This is why we assume in the following that the 
current j(x) is nonzero only inside the metal. 



Boundary conditions 

The magnetic field itself is continuous across the interface of 
the metal (assumed to be non-magnetic). The electric field 
components parallel to the interface are continuous as well. 
For the normal components, the continuity of the displacement 
field gives 

21 n Fl - (2?r)2 



A 2 



El 



(4) 



We shall assume that the wavelength A is much larger than any 
other relevant length scales in the problem, and take the limit 
5 <C A — * 00. To lowest order, this transforms Eq.(@) into 



the interface and find, using that the current j vanishes at the 
interface, 



V x B|° n ut = Mo aE|°f = 



(7) 



The last equality is again due to the negligible vacuum conduc- 
tivity. (Note that e^uja — (2ir6/\) 2 /2 — > in lowest order.) 
Taking components parallel to the interface, we find 



n x (V x B) 



-\i§a n x E 



(8) 



Now, from V • B = = 9jt(V • B) = did^Bi, we can derive 
using the Gauss theorem: rijVBi|™ = 0. This identity can- 
cels one of the terms coming from the expansion of the dou- 
ble vector product in 0. We thus find the jump condition 0. 
Note that the present derivation is valid for any geometry of the 
interface — which is less obvious for the previous one because 
one has to integrate the Laplacian operator. 



3 Incoherent summation 

In this section, we focus on a planar metallic surface and the 
quasi-static limit (geometrical distances even smaller than the 
skin depth). We display the surface electric field and the cor- 
rection it implies for the transmitted magnetic field. Technical 
details are deferred to AppendixlAl 



3.1 Transmitted field expansion 

For a planar surface parallel to the a;y-plane, an expansion in 
plane waves with two-dimensional wave vectors is straightfor- 
ward. We use the notation K = (k x , k y , 0) and find just below 
the metal surface the following expansion 



n-EL = 0. 



(5) 



B(r) 



d 2 K 



•BJK1 



(9) 



We check explicitly below that this boundary condition yields 
results consistent with a fully retarded calculation. 

We next calculate the jump condition for the magnetic field 
due to the surface current. Assume first a planar metallic sur- 
face located at z = with unit outward normal n. Noting that 
Vcr(r) = — and(z) and integrating the magnetic wave equa- 
tion along a path perpendicular to the interface, we find the 
following jump condition 



dn 



= [loan x E 



(6) 



where the scripts 'in' and 'out' mark the field inside and outside 
the metal. The electric field is taken at the interface because 
only its tangential components are involved. 

We can also get this boundary condition directly from the 
Maxwell equations and without specifying a planar boundary. 
At the metal surface, we evaluate just above and below 



where kj >r = K ± i/tn with k 



= yjK 2 - q 2 (Re k > 0). In 
the quasi-static limit, we have K ^> \q\, and therefore k « K. 
More general formulas can be found in AppendixlAl As shown 
there, a current density localized below the metal surface pro- 
duces an 'incident' magnetic field with Fourier transform 



B i [K] = ^k i xJ[K] 



(10) 



where 



,0 

J[K] = / dx 3 e Kxs j[K; X3 } (11) 
J — 00 



with j[K; £3] being the 2D spatial Fourier transform of j(x). 
For the solution of the reflection/transmission problem, we also 
need the electric field at the interface. Its tangential components 
are given by 

n x E[K] = - n * (kj • J[K]) =2nx E;[K]. (12) 
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The factor of two between the 'incident field' and the actual 
field at the interface, may be explained intuitively by working 
with image currents (or dipoles): to fulfill the boundary condi- 
tions, one combines the fields of the actual dipole (inside the 
metal) and of an image dipole. Since the field is incident from 
a region with a large 'refractive index', the reflected field has 
the opposite sign: source and image dipoles therefore have the 
same polarity if they are parallel to the interface. As a conse- 
quence, their field components parallel to the interface double. 
Dipoles perpendicular to the interface have mirror images with 
the opposite polarity, leading to the cancellation of the normal 
field component, as required by the boundary condition (0. 

Combining the jump condition (|6j for the magnetic field 
with the in-plane electric field dl2i . we get the transmitted mag- 
netic field: 

B t [K] = ^|k t x J[K] + ^5_(n x K)(k t • J[K]) (13) 

where k t = K + iA'n. With respect to the incident field dlOi . 
we thus have an additional term with components parallel to the 
interface. This term is absent when computing the field gener- 
ated by the current density j as if the latter were located in 
vacuum. It ensures in particular that the transmitted field van- 
ishes if j is parallel to n, which is a well-known result for a 
planar geometry (see, e.g., |4|). We show now that the interfer- 
ence between both terms reduces the field fluctuations for some 
components, and reproduces the noise tensor found asymptoti- 
cally from a retarded calculation. 

3.2 Magnetic field correlation tensor 

We first compute the correlation function of the Fourier trans- 
formed current. Local thermodynamic equilibrium gives the 
basic relation 1 13 23 1 for the current fluctuations 

Ofc(xi; w)jz(x 2 ;o/)) = 2%S(x.i;w)SuS(xi - x 2 )<5(u> - a/), 

(14) 

and we find 

{JJi\K 1 ;u]Ji\K 2 ;u']) = ^-S(u)8 k i8(K l -K 2 )S(u-u') 

(15) 

The noise spectrum behaves like S(x;ui) ~ 2a(x)kBT at 
low (sub-thermal) frequencies. The magnetic noise tensor 
Bij (xi , x 2 ) is defined by 

(B?(xi;o;)fl i (x2;w , )> = 1-k8(u> - u 1 ) B^xi, x 2 ) (16) 

Whenever no confusion is possible, we suppress the frequency 
dependence for simplicity. At low frequencies, the magnetic 
noise spectrum tends towards a constant anyway. 

The spatial Fourier transformed magnetic field Jl 3i thus 
yields, using the current correlation function \\5\ . 

B^xx.xaj 4 j 2{2n)2K5 e * v 

(17) 

Xij = (-\Ke lM k* k + (n x K)jfc*) 5 lq 

x (iKejpgkp + (n x K)jk q ) (18) 



This tensor can be worked out using the relations k x k* = 
(K + U^n) x (K - iKn) = 2L£Tn x K and |k| 2 = 2K 2 , valid 
in the quasi-static limit. We find 

X i:j = K 2 (2K 2 8 ij - hk* - 2(n x K),(n x K)j) (19) 

where the last term includes also the crossed correlations be- 
tween the two terms in Ea. (ll3> . leading to the minus sign. For 
the planar geometry, we can work out the integral over the az- 
imuthal angle. For simplicity, we focus on the noise tensor at 
the same position r = Xi = x 2 . The angular average (denoted 
by double brackets) gives 



(hk*)) = ^K 2 A tJ +K 2 n tnj 
«(n x K)j(n x K),)> = \k 2 A t] 



(20) 
(21) 
(22) 



where Ay = diag(l, 1,0) is the in-plane Kronecker symbol 
and Sij — diag(i,i,l) an anisotropic tensor that was also 
found inRef. OH, 

using the asymptotic expansion of the fully 
retarded magnetic noise tensor. Thanks to the additional term 
in ( I19> . we thus find the correct magnetic correlation tensor. 
Without this term, diag(|, |, 1) would have come out. Let us 
check the prefactor of the spectrum. It is given by 



Bij(r;w) 



til 
4 

tis_ 

32ttz' 



-2Kz 



' AK 

167TZ 



ix l3 )) 



(23) 



This is the result given in Eq.(24) of Ref. OH, taken in the 
quasi-static limit (distance z small compared to the skin depth). 
We thus have shown that when the correct boundary conditions 
for the magnetic field are used, the 'incoherent summation' ap- 
proach is equivalent to the more rigorous FD theorem. 

3.3 Impact of finite skin depth 

We now discuss what the previous formulas become when 
the relevant geometrical distances are comparable to the skin 
depth 6. The transmitted magnetic field takes the form (see Ap- 
pendixlAl 



BJK1 = 



Mo 



K 



k( (nx K) J[K 



(24) 



Here, K is the unit vector along K. This result is transverse as 
it should because the vacuum wave vector satisfies k 2 = 0. It 
also coincides with Ea.( ll3l in the quasi-static limit, as a simple 
calculation shows. Note that, again, B t = if j||n. 

The magnetic correlation tensor obtained from Eq.(I24> can 
be worked out and reduces to the following simple formula 



ti<rk B T f 00 K 3 dKe~ 2Kz 



2n 



Ren \k + K\ 



(25) 
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It is interesting to note that the anisotropy tensor Sy = 
diag(^, 5, 1) describes the magnetic noise through the entire 
distance range. The asymptotic limits for the integral are 



K 3 dKe~ 2Kz 

: ) Re K I K 



K 




(26) 



The second limit reproduces the result obtained asymptotically 
from the exact solution in the skin-dominated regime. 

Let us analyze this skin-dominated limit in more detail and 
consider here the case \q\ = 0(1/ 5) 3> K. We then have 
k w — iq 3> K and get to leading order 



B t [K] « ^k ( (nxK)-J[K] 



(27) 



Comparing to 1241 . we see that the skin effect effectively pre- 
vents the magnetic field from leaking out of the metal. It is 
sufficient to work to this order to get the large distance asymp- 
totics rEa. (l26> . second line]. 



3.4 Discussion 

We start with a discussion of the change in the power laws (I26> . 
as one changes the distance from below the skin depth to much 
larger values. In the short distance regime, the effective volume 
inside the metal that contributes to the magnetic noise, is of the 
order z 3 , since across the distance z, absorption is negligible. 
Adding incoherently magnetic fields with an amplitude ~ 1/z 2 
for each element in this volume, gives the 1/z power law for 
the magnetic noise power. At larger distances, damping in the 
metal becomes relevant, and one expects only a surface layer 
of volume ~ z 2 S to contribute. This leads to an scaling S/z 2 
that is not the one found here. In fact, as the skin depth gets 
shorter, the transmission through the metallic surface also be- 
comes more inefficient, as discussed in Sec 13. 31 This leads to a 
reduced transmitted field. The calculation shows that in Fourier 
space, one factor 1/K becomes 1/q [Eos. ( 1131 1271 1 so that the 
transmitted field scaling changes like 1/z 2 1— > 5/z 3 . We get 
the <5 3 /z 4 behaviour by incoherently summing this up over the 
near-surface volume ~ z 2 8. 

Let us compare to results obtained previously for metal- 
lic layers with finite thickness t. Calculations in this geometry 
have been performed by Varpula and Poutanen |4|, Sidles and 
co-workers |6|, and Rekdal and co-workers 1181 . At low fre- 
quencies where the skin depth becomes the largest scale, one 
gets 



5 > z.t 



16tt 1 



t 



z{z + t)' 



(28) 



This is consistent with the simple rule of removing the vacuum 
half space below the layer (replace 1/zby 1/z - l/(z + t)), 
ignoring the boundary conditions at the lower interface. One 
can actually show that sub-layer material add negligible noise 



as long as its conductivity is much less than that of the metallic 
layer (see Refs. 1181251 1. In the limit of a thin layer, t <C z, 
the noise is smaller compared to a half-space, because it de- 
cays like t/z 2 . For cylindrical wires, the noise reduction is even 
stronger, see Refs. 1171211 . 

At higher frequencies, the skin depth becomes shorter, and 
different regimes emerge depending on the relative magnitude 
of distance z and thickness t. We focus in the following on the 
thin layer limit. Varpula and Poutanen |4| give the following 
empirical interpolation formula for a layer thinner than the skin 
depth 



16tt 



t/z 2 



1 



[Azt/{-K 2 5 2 )} 2 

(29) 

This gives at large distance z^ 5 a. noise spectrum with a scal- 
ing ~ S 4 / (tz 4 ), similar to the half-space, but with an increased 
amplitude (by a factor of order 6/t). Note that in this regime, 
decreasing the layer thickness just produces the opposite effect 
on the noise. This unusual result has been confirmed exper- 
imentally in the kHz range |4|. A non-monotonic behaviour 
with either skin depth or conductivity a = 2/(fiQUjS 2 ) has 
also been pointed out by Rekdal and co-workers 11711 81 : bad 
conductors show only weak current fluctuations, while good 
conductors efficiently screen the magnetic field. 1 Eq.j29> in- 
deed yields a noise maximum for 5 ~ \[zt, consistent with 
Ref. IT8l . 

Finally, a layer thicker than the skin depth has been consid- 
ered in Refs. 1 18 6 1, where the following asymptotics is derived 




t < 6 < Vzi 



S <C min (z, t) 



(30) 

Note that the first line differs from Eq.j29> by a numerical fac- 
tor - this may be due to the chosen interpolation. The second 
line, consistent with Ref. 1 18 1, shows that for a very short skin 
depth, there is no difference between a metallic layer and a 
half-space rEa.J26H. This could have been expected given the 
highly efficient screening. 

Let us finally touch upon the case of a superconducting 
object. Sidles and co-workers 1 6 1 have argued that it suffices 
to use a complex conductivity and to make the replacement 
a 1 ► Re[<r(u;)]. For an ideal superconductor and zero temper- 
ature, London's equation yields <j(lo) — — iX\/(noUJ) with Xl 
the London penetration depth, and magnetic near field noise 
is completely suppressed. At finite temperature, the supercon- 
ducting phase coexists with a normal phase, and Re[cr(w)] is 
finite. In terms of the (frequency-dependent) phase angle ip in 
a = \a\c~ iip , the following interpolation formula is given in 



A similar situation occurs for the absorption of normally incident 
plane waves in a thin metallic film. The maximum absorption occurs 
for S ~ s/Xt. See 1 6 1 for the link to magnetic noise and, e.g., (24 1 for 
an instructive discussion. 
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Ref. |6| for the magnetic noise spectrum above a supercon- 
ducting layer 2 



V 



k B T | a | cos (p 36 3 t 
lfor SlJ ~D~ 
38 3 z(z +t) + 2^2(1 - e~ Qc ) x 
x tz 2 (z + SV2) 2 cos(vr/4 - ip/2) 

zt + AS 2 sin ip 
V2 z6 cos(ir/4 - tp/2)' 



(31) 

(32) 
(33) 



Numerically, it is found that this formula reproduces the results 
of an exact calculation to within 2 dB. For a normal conductor 
(<p = 0), it reproduces the asymptotics (1261 1301 . 

Sidles and co-workers have also given corrections for a 
material with a weak magnetic susceptibility — A*oj > 
0) where thermal magnetization fluctuations contribute to the 
magnetic field noise as well (with a noise spectrum propor- 
tional to Iml//i), see Eqs.(35, 36a) of Ref. |6|. Rekdal and 
co-workers |18| have pointed out that measurements of the 
magnetic susceptibility actually allow to infer the frequency- 
dependent complex conductivity that determines both the skin 
depth and magnetic noise properties. For niobium, the skin 
depth significantly differs from the London penetration depth at 
temperatures below the transition point and frequencies around 
500 kHz. 



4 Magnetic Green function 



In this section, we switch to an alternative approach to mag- 
netic near field noise that exploits a link to classical dipole radi- 
ation. In fact, the field radiated by a single point-like magnetic 
moment is sufficient to get the magnetic noise spectrum when 
the fluctuation-dissipation theorem J34I is used. This is a sig- 
nificant advantage for numerical calculations that are needed 
anyway in more complex geometries. In the 'incoherent sum- 
mation' technique, one not only faces a similar effort to be in- 
vested in the computation of the field, but the calculation has to 
be repeated for a large number of inequivalent sources (all vol- 
ume elements filled with absorbing material). As pointed out 
by Sidles et al. [6|, this redundancy can be avoided using the 
reciprocity theorem: once the field emitted by a suitable point 
source is computed, the relevant quantity is the total power ab- 
sorbed in the metallic structure. In the following, we formulate 
the equations to solve near arbitrary metallic structures, with 
retardation in vacuum being neglected. Subsection l4.2l special- 
izes to a half-space and shows that the reflected Green tensor is 
consistent with a fully retarded calculation. 



4.1 Fluctuation-dissipation theorem 

The fluctuation-dissipation (FD) theorem for the magnetic field 
reads 1 1 5 1 : 

B ii (x 1 ,x 2 ;o;) = 2hf(nu/k B T)lm7U j (xi,x.2-,w) (34) 

where H is the magnetic Green tensor, i.e., the magnetic field 
radiated by an oscillating magnetic point dipole m at X2, 
-Bdi p ,i(xi; uj) = J2j T~(-ij (xi, x 2 ; uj)mj. For the low-frequency 
regime relevant here, the temperature dependence reduces to 
f(hu/k B T) = k B T/hiu. 

Since we are interested in atoms trapped in vacuum above 
a metallic structure, we shall take xi = x 2 in vacuum. The 
magnetic dipole field Bdi P then can always be written as the 
sum of the vacuum radiation plus a field scattered or reflected 
from the structure. The vacuum field gives an imaginary part 
Im7i vac (xi, xi; u>) that reproduces Planck's formula for the 
blackbody radiation spectrum. We shall actually neglect this 
contribution compared to the one of the scattered field. For a 
planar surface, the scattered field can be written as an integral 
over Fresnel reflection coefficients (see, e.g., Ref. |16|). We 
check in the following section that the boundary conditions of 
Sec. |3 reproduce the Fresnel coefficients, at least in the low- 
frequency limit we focus on here. 

We shall use the vector potential A and the scalar potential 
<j> in the 'generalised Coulomb gauge' V ■ eA = 0. This gives 
the following wave equations for the domains outside and in- 
side the metallic objects whose shape is left arbitrary for the 
moment. Outside the object: 

V 2 <h = 



V 2 A = - 
V- A = 



HoV x M - — V0 



(35) 
(36) 
(37) 



where M(x) = m#(x — Xi) is the magnetization density for 
a point dipole. Inside the object: 



V 2 A 



V 2 c6 = 



V • A = 



(38) 

(39) 
(40) 



We have corrected an obvious error in Eq.(6a) of Ref. 1 6 1 and took 
the classical limit hui <C ksT. 



We will consistently work in the limit \qc/uj\ ~ A/<5 — > oo, 
with spatial derivatives being comparable to q. We thus cover 
length scales comparable with or smaller than the skin depth. 
Combined with the boundary conditions for the potentials, 
these equations allow to determine the field everywhere. Note 
that there is no source term in the equations for the scalar po- 
tential <fi, even when the boundary conditions are taken into 
account. Without loss of generality, we therefore put <fi = in 
the following. 

4.2 Planar geometry 

In the planar case, we have a simple analytical solution for the 
magnetic noise tensor - a benchmark result that has to be repro- 
duced by our theory. Details of the calculation can be found in 
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AppendixlBl We use the boundary condition n- A| in = char- 
acteristic for the metal-vacuum interface [Eq.Q] that follows 
from E = iojA. Translational symmetry allows to expand the 
field radiated by the magnetic dipole in vacuum in plane waves 
with wave vector K parallel to the interface. Focussing on an 
incident plane wave, the calculation yields a reflected magnetic 
field with Fourier amplitude 



B r [K] = §^k r (n \ K.) • A,[K]. 



(41) 



where A;[K] is given in Ea.( l53t . The ratio (K — k)/(K + k) 
is the same reflection coefficient that appears in the (retarded) 
magnetic Green function, when the limit A — > oo is taken, see, 
e.g., 1 16 1. In AppendixlBl the magnetic Green tensor computed 
from J41t is found to be 



B ii (r;w) = - 



^ 

2t:lu 



• K 3 dK2Imne- 2Kz 
\k + K\ 2 



(42) 



This expression agrees with Ea. J25t thanks to the identity 

2 Im k Re k = Im n 2 = — Im q 2 = — a aaj. 



5 Summary and conclusion 

We have discussed in this paper calculations for low-frequency 
magnetic noise fields at sub-wavelength distances to metallic 
objects. The role of the object surfaces has been clarified: they 
screen some field components so that only current elements 
parallel to the surface produce fields outside the object. This 
occurs even on a distance scale where dissipation in the metal 
is negligible. Neglecting this effect leads to errors up to a fac- 
tor of three for the components of the magnetic noise tensor at 
short distance (smaller than the skin depth S), and completely 
wrong power laws at larger distances (3> 5). As a consequence, 
the simple incoherent addition of thermal noise fields has to be 
replaced by a more involved calculation, preferably based on 
the fluctuation-dissipation theorem for the field. We have for- 
mulated an outline of this calculation in a generic geometry, 
spelling out the boundary conditions that apply in the low- 
frequency regime characteristic for miniaturized atom traps. 
We hope that this opens a way to accurate and numerically ef- 
ficient methods of characterizing magnetic noise spectra near 
complex atom chip structures. 

Our theory can be extended to dielectric objects as well. If 
absorption is large and |e| ^> 1, the same approach can be car- 
ried over, with the skin depth defined by 1/5 = (u>/c) Im y/e. 
For a purely real permittivity, however, one has to include re- 
tardation in the vacuum regions to get a nonzero imaginary part 
in the magnetic Green function. When e is of order unity, the 
boundary conditions for the fields assume, of course, their stan- 
dard form. Numerical calculations are currently under way to 
test the validity of the non-retarded approach. 



This work was supported by the European Commission in the project 
ACQP (contract IST-2001-38863) and the network FASTNet (contract 
HPRN-CT-2002-00304) . 



A Transmission through a planar interface 

We first compute the magnetic field radiated by a current dis- 
tribution j(x) localized inside a homogeneous metal. This is 
given by 



JW 



(43) 



where we recall that q 2 = 2i/<5 2 . To solve the transmission 
problem through a planar interface, it is expedient to use the 
expansion of the spherical wave in plane waves (the Weyl an- 
gular spectrum) with the wave vector K = (k x , k y ,0): 



Jq\r- 



4tt r - 



d 2 K 
(2 7 r) 2 K < 



,iK-(r— x) — k\z — xa\ 



(44) 



If we consider a point r just below the interface, the absolute 
value in Ea.J44> can be dropped, and we get a Fourier coeffi- 
cient 

iMo, 



Bi[K] - ^k, x J[K] 



(45) 



where k s ; = K + i/tn, and J[K] is given by Eg. (II 1>. with K 
in the exponential replaced by k. Only in the quasi-static limit, 
K — > K, and we recover Eq.dlOli. 

In the same way, we get the normal derivative 



— Bi[K] = -kB s: [K] = -^k s: X J[K] 



(46) 



The first equation reflects the property that the magnetic field 
created by the current is propagating (in fact, decaying) in the 
upward direction. For the reflected and transmitted fields, sim- 
ilar relations hold: 



d 

— B r [K 

On 







-AtB r [K]j — B t [K 
on 



-KB t [K\. (47) 

The jump condition l|8j thus gives 

- KB t [K] + k (Bj[K] - B r [K]) = ^ra x E t [K]. (48) 

It turns out that to proceed, we do not actually need to com- 
pute the electric field at the metallic surface. We combine the 
continuity of the magnetic field perpendicular to the surface 
with the corresponding component of Ea.J48> and solve for the 
transmitted field component 



n B t [K] = 



2 k 



K 



n • Bj[K] 



(49) 



I thank Isabelle Bouchoule for stimulating comments and the Labo- 
ratoire Charles Fabry of the Institut d'Optique for its kind hospitality. 



where we recognize one of the Fresnel coefficients. The com- 
ponent along K follows from div B t = ik t ■ B t [K] = 0. 
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The last component is along n x K, and we can use the 
following trick to show that it actually vanishes. The Maxwell 
equations yield (a = i,r,t) 

i(nxK)-B Q [K] = in-(k Q xB Q [K]) = ^ CT Q n-E Q [K] (50) 

with (Ti >r = a and at — > 0. We now use the boundary condition 
that the normal electric field is zero inside the metal [Eq.(|5j], 
and get from Eq.(|50j (n x K) ■ (B 4 [K] + B r [K]) = 0. Now, 
this is a magnetic field component tangential to the surface, and 
therefore 

(nxK).B f [K] = 0. (51) 

Combining the two nonzero components found above, we read- 
ily get the transmitted field (1241 . 



B Magnetic Green tensor 

The vector potential created by a magnetic point dipole in free 
space (solution of Ea. (l36» is of the form 



A i( r) = V r x 

47r r — xi 



(52) 



where xi is the dipole position. Above a planar surface, we use 
the Weyl expansion J44l > and find the following Fourier coeffi- 
cient for the field incident on the interface: 



AJK1 = -^k 4 x me- ik *- Xl 
2K 



(53) 



with an incident wavevector k, = K — iK n. The reflected field 
is characterized by a wavevector k r = K+LftTn and transversal 
as well, i.e. k r • A r [K] = 0. 

The transmitted field has the wavevector k f = K— i«n with 
K 2 — k 2 = q 2 (Re n > 0). Here, we have two 'transversality 
conditions': 



k t -A t [K] = 0, n-A t [K]=0, 



(54) 



where the second one actually comes from the boundary condi- 
tion for the electric field, Eq.(|5|i. We conclude that the transmit- 
ted field is parallel to the vector n x k t = n x K. This vector 
lies inside the boundary and is perpendicular to K. Since the 
tangential components of the vector potential are continuous, 
the reflected field A r has to cancel the component of A; paral- 
lel to K. This gives a first condition for the reflected field: 



K-(Ai[K] + A r [K]) = 0. 



(55) 



We need a second boundary condition to solve the problem. 
This is, of course, the continuity of the magnetic field B = 
V x A. The magnetic field on the inner side of the interface is 
computed to be 



B 



ik t x A t = i (nK 2 + ikK) t 



(56) 



where t is an un-normalized transmission coefficient. Comput- 
ing the normal and K component on the 'outer side', we get 
the linear system 



K 2 t = -in • B| out - (n x K) • (A, + A r ) 



(57) 



iK B| 



iK (n x K) ■ (A; - A r ) (58) 



whose solution involve the standard Fresnel reflection and 
transmission coefficients: 



t = 



(n x K) • A r 



K{K-\ 

K-k 



(n x K) • Ai 



(n x K) • A, 



K + K 

This yields the following expression for the reflected field 

K 



(59) 
(60) 



K 



(K • A,) (n x K) (n x k) A, (61) 



whose first ('longitudinal') term ensures condition J55i while 
still being 'transversal' (this is due to the fact that k^ = in the 
limit A — > oo). The corresponding magnetic field is determined 
by the second term only and one finds Ea.J41>. 

From Eq.( I61> . we identify the following magnetic Green 
tensor 



H{r, r) = H (vac) (r, r) + H (rof) (r, r) 
d 2 K Li e~ 2Kz K- k. 



W (ref) (r,r) = 



(27T) 



2K K 



(62) 
(63) 



using the identity (n x K) x k^ = ik^. The integration over the 
azimuthal angle amounts to angular averaging: ((k r <£> k«)) = 
K 2 Sij. For the fluctuation-dissipation theorem, we need the 
imaginary part of the reflection coefficient 

K-k —2K Im k 



Im • 



K 



(64) 



The vacuum Green tensor is purely real in the static limit A — > 
oo we focus on here, and does not contribute to the magnetic 
noise spectrum J34J> . Putting everything together, we get the 
Green tensor J42I . 
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